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ABSTRACT 

In this Letter we present calibrated observations of the average mass of galax- 
ies with 22 < / < 25 in the Hubble Deep Field. We measure the mean mass 
profile using the statistical gravitational lens distortion of faint blue background 
galaxies, based on 2221 foreground-background pairs. The observed weak lens- 
ing distortion is calibrated via full 3-D simulations and other HST data on a 
foreground cluster of known mass. 

Inside a projected radius of 5" we find a 3cr detection of shear. Fitting to 
an isothermal mass distribution, we find a shear of 0.0661oq]J at a radius of 2" 
(~ 8 h^ 1 kpc). The (la) limits on the parameters for the truncated isothermal 
model are a v = 18513° km s _1 , r outer > 15 h^ 1 kpc, and r core = 0.611q I hT 1 kpc. 
This corresponds to an average galaxy mass interior to 20 h~ l kpc of 5.9j|;7 x 10 11 
h^ 1 Mq. Inside 10 h^ 1 kpc, we find an average rest-frame V band mass-to-light 
ratio of 11.4±|;? h (M/L V ) Q . 



Subject headings: dark matter - galaxies: fundamental parameters - gravitational 
lensing 



1. Introduction 

The mass in the dark halos around individual field galaxies has long been debated. While 
there are compelling kinematic data on our own galaxy (Fich & Tremaine 1991) and pairs 
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of galaxies (Zaritsky and White 1994), suggesting halos with equivalent circular velocities of 
~200 km sec -1 extending as far as 30 kpc, direct measures of mass using weak gravitational 
lens distortions have been difficult. Typical galaxies produce ~ 1" deflections, resulting in 
relatively small shear at radii > 5". Strong lensing of background QSOs has provided a 
spot check of total enclosed mass in several cases (Kochanek 1991), confirming masses of 
~ 1 — 2 x 10 11 h~ l M within 4 h~ l kpc. However, the known lensed QSOs may selectively 
sample the tail of the galaxy mass distribution (Wallington & Narayan 1993). 

Statistical gravitational lens distortion in a large sample of galaxies can address the 
more general question of the mean mass of all galaxies. Previous attempts to measure the 
average total mass of a large sample galaxies by statistical weak lensing were hampered by 
seeing and systematic error (Tyson, et al. 1984 [TVJM84]; Tyson, 1987). Studying weak 
shear at radii > 5" in a 4.8' radius field taken with the Palomar 5-meter telescope, Brainerd 
et al. (1996) recently claimed evidence for massive halos extending beyond 30 hr l kpc. The 
Hubble Deep Field (HDF) data, being multi-band, high spatial resolution, and very deep - 
despite the small size of the field - provide a large volume sample and an opportunity to 
control systematics and bring the <5" mass signal out of the noise. 

Here we report a calibrated measurement of the mass profile of field galaxies in the 
HDF. Galaxies with the faint surface brightness and colors of the distant faint blue galaxy 
population are used as background galaxies. 2221 foreground-background pairs of galaxies 
are examined for the induced shear from the foreground galaxy. We compare our result for 
the HDF with full 3-D simulations, with deep HST weak lens studies around a cluster of 
galaxies, and with broader but "shallower" HST surveys. A comprehensive description of 
the cluster HST WFC2 weak lens shear efficiency calibrations, and comparison with other 
galaxy-galaxy weak lensing observations using HST and ground-based imaging will be given 
elsewhere. 



2. Image Processing and Photometry Catalogs 

For this study, we used the "drizzled" HDF version 2 data in which the dithered images 
were combined on a resampling scale of 0.04" /pixel. The HDF observations and pipeline 
processing are described in Williams et al. 1996. These reconstructed images cover 4.3 
arcmin 2 for all three WFC2 chips. The sky noise on 1" scales in these reconstructed images 
is 28 B mag arcsec -2 and 28.5 R mag arcsec -2 . For this analysis we limit ourselves to the three 
deepest bands: F814W, F606W, and F450W (hereafter called I, R, and B). We re-analyze 
the images using FOCAS to produce photometrically accurate matched catalogs. The zero- 
points in the various filter bands are set by the HDF calibration (Williams et al. 1996). 
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We calculate magnitudes and colors using an equal-isophote method: initial object 
lists are constructed by running FOCAS on B + R + I master images. This produces 
a master list of objects and isophotal regions. Magnitudes are then computed from the 
individual filter images using these identical isophotal areas, so the color estimates are not 
biased by systematic variations of the isophotes from band to band. This procedure yields 
magnitudes and intensity-weighted moments measured accurately in three bands. At a 
limiting magnitude of 29 B mag, about 40 galaxies between 24 and 29 mag will lie within a 
10" circle projected on any foreground galaxy. 

In order to extract the lensing shear, we must reliably measure the ellipticity and ori- 
entation of a faint galaxy. The faint background galaxies we use are quite small: the profile 
of a typical 27th mag arclet is shown in Figure 1. Useful surface photometry of these faint 
blue galaxies in the HDF data is possible only within r ~ 0.5 arcsec. Simulations of random 
galaxy shapes using the observed noise characteristics of the HDF field show that we can 
measure ellipticity for galaxies as faint as 28 B mag and ellipticity 0.3 with an uncertainty 
in ellipticity less than 0.1. This suggests that we can safely use galaxies down to B = 27.5 
as potential probes of the shear. 



3. Arclet Analysis and Mass Model Simulations 

To minimize the contamination from foreground galaxies, we use the color and appar- 
ent magnitude of the galaxies as rough statistical redshift estimators. Any galaxy in the 
magnitude range 22 < / < 25 (roughly redshift 0.2 - 0.7) is a candidate foreground lensing 
galaxy. Blue galaxies with B — R < 0.3 that are fainter than a candidate foreground galaxy 
but brighter than 27.5 B mag are candidate background galaxies. Due to competing volume 
and faint-end luminosity function effects, these faint galaxies cover a wide range in redshifts 
extending from z = 1-3; blue low surface brightness sources are likely to be at the highest 
redshifts. The trend towards redshift 1 at 25th mag is clear in a magnitude - \og(z) plot 
(Koo & Kron 1992; Schade et al. 1995; Lilly et al. 1995). 

Using this selection, we find 110 candidate lenses and 645 candidate background objects 
yielding 2221 candidate foreground-background pairs in the three WFC2 CCD area of the 
HDF. We use FOCAS to calculate the intensity-weighted second moment for each source 
galaxy image (arclet) within 15 arcsec of a candidate foreground galaxy. We calculate the x 2 , 
y 2 and xy intensity moments (apodized with a Gaussian kernel of size equal to the half-light 
radius) transformed to (r, 9) moments relative to each foreground galaxy center (TVJM84; 
Tyson, Valdes and Wenk 1990). These principal- axis transformed second moments of the 
arclets within some annulus about all candidate foreground galaxies are averaged together 
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to give the observed shear at that radius. 

This observed shear as a function of radius for all 2221 pairs is plotted in Figure 2. 
The points are the observed independent shear values in each radius bin, and the error bars 
obtained via the full simulations are consistent with bootstrap resampling errors. As always, 
the largest noise source is the randomly oriented non-zero ellipticities of the source galaxies. 
We find a 3c significance positive mass effect observed at radii less than a few arcseconds 
(~ 10 h~ x kpc for 0.2 < Zf g < 0.7). While we have not "tuned" the arclet color to maximize 
shear, if we relax the blue selection this shear decreases; this is consistent with the blue B-R 
colors of the highest redshift galaxies. 

What galaxy mass distribution would generate this observed shear? In selecting our 
sample galaxies based on brightness and color, the observed shear will be diluted by the 
convolution of the separate redshift distributions of the candidate lenses and arclets. The 
WFC2 point spread function (PSF) and the drizzle algorithm have an effect on observed 
shear for arclets whose surface brightness falls below about 28 mag arcsec -2 inside 0.2" 
radius. Therefore we must calibrate the lens mass which would generate this observed 
statistical shear. We do this in two steps: (1) full 3-D simulations of the HDF, and then (2) 
comparison of a similar 3-D simulation with the weak lensing shear of this same faint blue 
galaxy population observed with HST WFC2 around a cluster of galaxies of known redshift 
and mass. 



3.1. Mass calibration 

We calibrate the shear-mass relationship by simulating galaxy-galaxy lens distortions 
including multiple scattering effects. Realistic simulated field galaxies consistent with mild 
luminosity evolution (Im et al. 1995) are randomly distributed in seven redshift shells. These 
simulated galaxies have the same properties as the observed galaxies (number- magnitude, 
size-magnitude, ellipticity, and redshift distributions) as found in recent deep redshift and 
angular size surveys (Schade et al. 1995; Mutz et al. 1994). The redshift distribution we use 
is similar to that derived from four-color photometry of the HDF (Mobasher et al. 1996), 
although lacking a distinct peak at z ~ 2. 

The lens distortions for galaxies in each redshift shell are computed using the galaxies 
in all of the shells in front of it. After the cumulative distortions for all redshift shells are 
calculated, the image is convolved with the PSF, binned down to the CCD resolution, and 
HDF-derived photon and readout noises are added. All galaxies are modelled as soft core 
isothermal distributions (Grossman & Narayan 1988), but with an outer mass cutoff, giving 
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for the spherical case a surface density distribution outside the core: 

S = S (/3/2r)(l + r 2 /rL er ) _1 (1) 

where E = a 2 /2G, a v is the line-of-sight velocity dispersion, r outer is the outer mass cutoff, 
and (3 = (1 + r^ ore / rl uter ) . For the candidate background galaxies, we measure the image 
distortion in terms of the (r, 6) second moments i rr and igg: 

T , v = iee - irr = l(r) 
i ee + i rr 1 - k(t) 

where the convergence n(r) = S(r)/E c and the shear 7(7") = [E(r) — E(r)]/E c (Miralda- 
Escude 1991), and where E c is the critical surface mass density. Outside the core, k << 
1, and we shall call T(r) the "shear". We plot curves for the resulting simulated shear 
T(r), corrected for the cluster calibration and efficiency described below, for three values of 
foreground lens galaxy mass in Figure 2. All three have outer cutoff radii of 30 hr x kpc, soft 
core radii of \XShr l kpc, and line-of-sight velocity dispersions as labeled. 

The critical surface density is related to the distance ratio: E c = c 2 / (4nGD). The 
distance ratio for a foreground-background pair is (cf. Blandford and Narayan 1992) 

_ (1 - q - did^jdi - d 2 ) 
~ (l-q -d 2 )(l-d 2 )(l + z fg ) 

where d 1 = y/1 + q Zf g and d 2 = y/l + q z bg- By averaging over many different galaxies 
we are averaging D over the distributions of Zf g and zj, g . We combine these uncertainties, 
including multiple scattering and the FOCAS arclet detection efficiency, in a single "efficiency 
function" Q(r). The value of Q(r) is determined by comparing the prediction for single 
galaxy-galaxy lensing from equation (2) with the results of the 3-D simulations of the HDF 
galaxy-galaxy lensing, including multiple scattering, normalizing to the galaxy cluster data. 
In this way, we construct an absolute calibration for the shear-o^ relation. A total of 4368 
CPU-hours of time on an SGI Power Challenger were needed to perform the simulations. 
Figure 3 shows a plot of Q versus r. Over the radial range of our data a constant Q = 0.33 
is a good overall approximation. 

If we have used the correct redshift distributions, etc., in our simulations, this procedure 
should give reliable shear predictions. We can largely eliminate such redshift uncertainties 
by simulating a lens of known mass. To fix the mass zero point for the simulation-based 
mass calibration we cross calibrate to a deep HST image of a massive cluster of galaxies at 
z = 0.4. The same 3-D simulation is run for CL0024+1654 and compared with the observed 
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weak shear in the HST WFC2 image of this cluster, using identical arclet selection rules. 
We find a similar (but 10 percent lower) efficiency for shear measurement in the cluster data 
at a radius where the mean shear is equal to the HDF shear at 2". The curves plotted in 
Figure 2 have these calibration corrections applied. 

3.2. Systematics 

A minor correction is made for the field distortion and the spatial variation of the 
WFC2 PSF by using star-dominated WFC2 images of the LMC. Corrections due to the 
drizzle algorithm were found to be quite small: measurement of the signal using undrizzled 
data yielded statistically indistinguishable results. Other systematic errors which would be 
present in the HDF data, via photometry or sample selection or shear analysis, are present 
also in the full simulations. Therefore, by comparing observed shear in the PSF-corrected 
data and simulations, the derived limits on lens mass are not biased by systematics. One 
sanity check is to reverse the pairs: use arclets as lenses and lenses as arclets. This gives 
a zero shear (0.006 ± 0.019). Another check is to scramble the coordinates of arclets in 
the catalog and recalculate the shear: again null (0.009 ± 0.02). Similarly, scrambling the 
positions of the foreground galaxies yields a zero signal (0.01 ± 0.02). Finally, each of the 
three WFC2 chips independently shows positive shear within 2" in the HDF. As a check on 
the FOCAS software, we measured the shear signal from a set of simulations identical to the 
ones used for calibration but with the lensing turned off. This resulted in a null signal (1.2a 
negative in the innermost bin). 

However, spiral arms in foreground galaxies in the HDF might occasionally be split off by 
FOCAS and mistaken for arcs. The simulations include the ellipticity distribution of bright 
galaxies, but not spiral arms. To check this, we take wider field deep ground-based imaging 
(0.5" pixels) and run FOCAS with the same parameters as with the HDF data, calling each 
pixel 0.04" and scaling the magnitudes accordingly. This generates 38900 pairs of galaxies 
with scaled apparent magnitudes and angular sizes similar to the HDF. We obtain a null 
result for the observed shear in this surrogate data: at 2" radius (25" in the original wide- 
field data) we find a mean shear of 0.0011 ± 0.0016, and at 1" we find 0.0019 ± 0.0025. Thus, 
spiral arms are apparently not a source of systematic error in the HDF shear measurements, 
down to 1" radius, on the scale of our observed shear. 
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4. Average Galaxy Mass 

Using the product of efficiency and average distance ratio Q(r) found in the cluster- 
calibrated 3-D simulations to normalize the analytic expression in equation f to the HDF 
shear data, we find a best fit galaxy velocity dispersion a v of 185^3° km sec" 1 for outer 
truncation radii larger than 15 h~ x kpc and r core < 2 h~ l kpc. This corresponds to a shear 
at 2" of 0.066l8;8iE- The average mass interior to 20 h' 1 kpc is 5.9l|;f x 10 11 h~ l M . For 
a v = 185 km sec -1 , the core radius is r core = 0.6±qJ kpc. From the absence of an odd 
image in strongly lensed QSOs, Wallington & Narayan (1993) argue for mass cores < 200 
pc. The HDF data are consistent with this requirement, although our upper limits are not 
as stringent. 

Fixing a v at 185 km sec -1 , we find an acceptable fit for truncation radii r outer > 3.5" 
(> 15 hr 1 kpc). Inside rs 10 hr 1 kpc, the rest-frame V band mass-to-light ratio is 11.4l 6 ; 7 
h (M/Lv)q- VLA-based lensed QSO data sample smaller galaxy radii and find correspond- 
ingly smaller M/L (see Fassnacht et al. 1996). Our M/L for the whole field galaxy population 
suggests a halo increasingly dominated by dark matter beyond 10 h~ l kpc. In this Letter we 
have only examined isothermal models. The radial profiles of our observed and calculated 
shear are different; we will explore this at greater length elsewhere, combining various deep 
HST and ground-based data. 

We briefly compare our results with the recent observations of galaxy-galaxy lensing by 
Brainerd et al. (1996). Although both studies probe galaxy-galaxy lensing, they do so on 
different scales. The Brainerd et al. observations sample a region outside 5" or ~ 15 h~ x kpc 
for their sample. By contrast, our sample shows strong shear precisely inside this region. 
Extrapolating their shear at 10" to smaller radii yields a somewhat greater shear than is 
seen in the HDF. Ultimately, the distribution of galaxy mass on these larger scales will be 
addressed by comparably deep surveys covering hundreds of times the area of the HDF. 
In particular, these surveys will have to go significantly deeper than the "Groth Strip" 
which covers 144 square arcmin, imaged with the HST in two filters: F555W and F814W. 
This covers 27 times the area of the HDF, but 3 magnitudes shallower, resulting in fewer 
close pairs. Because the candidate arclets are brighter and hence lower redshift, the lensing 
efficiency is also reduced. For the Groth strip, we find a shear 0.015±0.009 for 492 foreground- 
background pairs within 5 arcseconds, roughly consistent with the HDF results, although 
with much greater uncertainty. 

We gratefully acknowledge the help of Tom Duff, Greg Kochanski and Rick Wenk, 
discussions with Andy Fruchter, Allen Mills, and Penny Sackett, and suggestions of the 
referee. Phil Fischer kindly supplied the WFC2 LMC data. We thank Bob Williams and 
the HDF team at STScI for rapidly making these data available to the community. 
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Fig. 1. — Surface brightness profile of a typical 27 mag arclet in the HDF. Accurate second 
moments for shear measurement require a limiting surface brightness of 28-29 mag arcsec -2 
for these galaxies. Note that because the surface brightness is measured by averaging in 
elliptical annuli, adjacent annuli sample some of the same pixels, and hence the data points 
are correlated and the error bars only represent the photon statistics. The dashed line is the 
average R-band WFPC2 PSF. 
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Fig. 2. — The observed shear of the background galaxies, in 1.5" independent bins, around 
the positions of foreground galaxies for 2221 foreground-background pair candidates in the 
three WFC fields of the HDF. The background/foreground selection was based on colors, 
surface brightness, and total magnitude difference. The lines show shear from full 3-D 
calibrated simulations of truncated isothermal models with a v as labeled, r outer = 30 h~ x 
kpc, and r core = 1 h^ 1 kpc. 
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Fig. 3. — The product of distance ratio and efficiency for recovery of shear as a function of 
radius from the foreground galaxy, from the simulations. Q includes the effects of source 
and lens redshifts, and is calibrated by normalizing Q S im/Qobs to a lens of known mass, 
CL0024+1654. The la error bars are derived from the scatter in measured shear in the 
simulations. The dashed line is the value of Q used in plotting the curves in Figure 2. 



